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Abstract 

We discuss the symmetry decomposition of the average density of states for 
the two dimensional potential V = x 2 y 2 and its three dimensional generalisa- 
tion V = x 2 y 2 + y 2 z 2 + z 2 x 2 . In both problems, the energetically accessible 
phase space is non-compact due to the existence of infinite channels along the 
axes. It is known that in two dimensions the phase space volume is infinite 
in these channels thus yielding non-standard forms for the average density of 
states. Here we show that the channels also result in the symmetry decom- 
position having a much stronger effect than in potentials without channels, 
leading to terms which are essentially leading order. We verify these results 
numerically and also observe a peculiar numerical effect which we associate 
with the channels. In three dimensions, the volume of phase space is finite 
and the symmetry decomposition follows more closely that for generic poten- 
tials — however there are still non-generic effects related to some of the group 
elements. 

PACS numbers: 03.65.Sq, 11.15.Kc, 03.50-z, 11.30.Ly 
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INTRODUCTION 



The role of chaos in the classical Yang-Mills fields has been examined by several authors, 
the studies typically being divided into two regimes. In the first, one studies the full field 
theory |]J and tries to determine such global measures of chaos as the spectrum of Lyapunov 
exponents and spatial-temporal correlations. In the second, one studies the homogeneous 
or zero-dimensional limit of the problem M which admits a more microscopic analysis. 
Following this approach, one is led to consider the three dimensional potential V = x 2 y 2 + 
y 2 z 2 + z 2 x 2 and its simpler two dimensional cousin V = x 2 y 2 . The two dimensional problem 
has been independently studied since it is an interesting dynamical system in its own right. 
Until Dahlqvist and Russberg showed otherwise f3j, it was commonly believed that the 
classical motion in this potential was completely chaotic. Although this is not true, it 
remains one of the most chaotic potential systems known. It also serves as a useful example 
of intermittency f|||. Far from the origin, the motion is confined within one of four channels 
within which the problem is adiabatic so that a trajectory behaves in a smooth, regular 
manner. Upon exiting the channel, the trajectory undergoes a burst of strongly irregular 
motion before re-entering one of the channels. This form of dynamics, regular behaviour 
with episodes of irregularity, is called intermittency and is found in various physical systems 
including the classical helium atom |J] and the hydrogen atom in a strong magnetic field 
0. The first of these is governed by a potential very similar in form to x 2 y 2 ||. The three 
dimensional problem shares the properties of strong chaos and intermittency although this 
has been less extensively studied. 

We will be interested in the requantisation of these potentials, particularly in their den- 
sities of states. As proved by Simon || and later analysed in greater detail by Tomsovic 
T0|| ) t ne t wo dimensional potential has a discrete quantum spectrum in spite of having an 



energetically accessible phase space of infinite volume. This potential therefore violates the 
semiclassical relation that the average number of quantum energy levels below energy E is 
proportional to the volume of energetically accessible classical phase space. In this paper 
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we discuss a related property of this potential - the manner in which the average density 
of states decomposes among the various irreducible representations (irreps). Normally, the 
ratio of the number of states belonging to a given irrep R of dimension cIr is roughly d R /\G\ 
n| , [i2|| , where \G\ is the order of the group. There are then small h corrections depending on 



the symmetry properties of the irreps JL2|| . We will show here that for the potentials men- 
tioned above, the symmetry "corrections" can be anomalously large and in two dimensions 
are essentially leading order in their effect. 

The relevant symmetry groups for the two and three dimensional potentials are 
and the extended octahedral group respectively ( "extended" because we allow for inversions 
as well as rotations.) These groups have 5 and 10 conjugacy classes of group elements, 
respectively, and we need to analyse them all in order to calculate the average density of 
each irrep. The method for doing this when there are no channels was discussed in Ref. Jl3| 



for reflection operations and JT4J in the context of the permutation group of symmetric 



groups. It was then developed in a more general context in Ref. ||12|| . For some of the classes 
which appear here, the analysis is a straight-forward application of this theory. For other 
classes, however, the channel effects make it inapplicable and we use a different analysis 
based on the adiabatic nature of the Hamiltonian deep in the channels as introduced in 
HI . In both two and three dimensions, each channel calculation involves an analysis of the 



subgroup which leaves that channel invariant. 

The structure of the paper is as follows. In the next section, we review the formalism 
used in constructing the average density of states from approximations of the heat kernels. 
The approximations are based on Wigner transforms of the Hamiltonian and of unitary 
transformations which correspond to the group elements. This formalism will be used in 
the central region of the potential but will be adapted for application to the channels. In 
section II we apply this to the two dimensional potential and show that there are very strong 
effects arising from this decomposition - much stronger than what one would expect for a 
normal bound potential. In section III, we verify these results numerically and also point 
out the existence of a subtle numerical effect which is only apparent on doing the symmetry 



decomposition. In section IV we introduce the three dimensional generalisation and discuss 
the Wigner transforms corresponding to the various group elements. In section V we do 
the channel analysis of the three dimensional problem and use this to get the final results 
for all classes. In three dimensions, the channel effects are less dramatic but still introduce 
modifications to what one expects for a generic potential. 



I. FORMALISM 

We will interest ourselves in the smooth average part p(E) of the density of states, often 
called the Thomas- Fermi term. There is also an oscillating part posc(E) given by periodic 



orbits |L5[ but we will not discuss this in great detail so in what follows we suppress the bar 
on the smooth functions. The specification of only concerning ourselves with the Thomas- 
Fermi term in the density of states is made by invoking h expansions rather than expansions 
involving oscillatory functions of 1/h. One way to find the Thomas- Fermi density of states 
is to work with the partition function (often called the heat kernel), which is the Laplace 
transform of the density of states, 

Z(P) = Tr(e-^)=C(p(E)). (1.1) 

In the presence of a symmetry group, each quantum state will belong to one specific ir- 
reducible representation of that group so we will consider the heat kernels of each irrep 
separately, 

Z R {0) = Tr (P R e-^) . (1.2) 



Pr is the projection operator onto the irrep R and for a discrete group is given by JTB 

PR = ^iT,Xn(9M9)- (1-3) 

M g 

The sum is over the elements of the group, \G\ in number, Xr{9) is t ne character of group 
element g in irrep R, cIr is the dimension of irrep R and U(g) is the unitary operator 
corresponding to the element g, 



(r\U(gm = (g- 1 T\iP) = iP(g~ 1 T). (1.4) 



One standard way to proceed [17] is to find the Wigner transform of the operators e 1311 



and Pr and integrate them to evaluate the trace. The Wigner transform Aw {<l,p) of a 
quantum operator A is a representation of it in classical phase space and is defined as 



-V(q,p) = J dx(q_+ | 
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q _±) e -Wa (L5) 



To leading order in ft, it is valid to replace ( e ~^) w/ by e~^ Hw where the Wigner transform 
of the quantum Hamiltonian is just the corresponding classical Hamiltonian. Traces are 
simply evaluated in this representation since 

Tr(i) = J dqdpA w (q,p) 

Tr(AB) = — ^ J dqdpA w (q, p)B w (q, p), (1.6) 

where n is the dimension of the system. As we will see below, naive application of these 
formulas may diverge in the channels; nevertheless, the formalism can be adapted. 

In the evaluation of the Wigner transform of the projection operators (|1.3|), we need the 



Wigner transforms of the unitary operators U (g) . This is discussed in detail in Ref. |L2 
the results for all possible group elements in two dimensions are, 

(f>(/)) H/ (q,p) = l 
(f/(crj)) w (q,p) = Tih5(q i )5(p i ) 

(U(Rg)) (q,p) » -^-Mqi)S(q 2 )5( Pl )5(p 2 ). (1.7) 
sin \ 2 ) 

The Wigner transform of the identity operator gives unity; the transform of a reflection 
operator gives the delta functions of the position and momentum corresponding to the 
symmetry plane; and, the transform of a rotation gives the delta functions evaluated at 
the symmetry axis. (The third result is exact for 9 = it, otherwise it has higher order h 
contributions.) An additional useful property is that the Wigner transform of the product 
of two commuting operators is simply the product of their respective transforms. Using this, 



we obtain from ( |1.7|) the following relations for the three dimensional operators 



(U(a 1 a 2 a 3 )) w (q, p) = 7T 3 h 3 5(q)6(p) 



U{oRe)) w (q, P) « ^-5(q)5(p). (1.8) 
sin 

The first of these says that the transform of the product of three perpendicular reflections 
gives delta functions in all coordinates and momenta. The second says that the transform 
of a reflection through a plane times a rotation about the perpendicular axis gives the same 
delta functions. In both (|1 . 7|) and (|1.8|) , the relative power of h equals the co-dimension of 
the set of points left invariant by the group element. We follow Ref. jl2| in constructing 
"class heat kernels" 

Z(g;i3)=Tr(u(g)e-^) (1.9) 

so that 

z R ([3) = ^Y.Xn(g)z(g;{3). (1.10) 

M g 

The functions defined in Eq. Q1.9Q are "class functions"; they do not depend explicitly on 
the group element g but only on the class to which it belongs. 



II. THE POTENTIAL V = x 2 y 2 

The equipotential curves of this potential are shown as the light curves in Fig. [l|. The 
symmetry group is C± v , the same as that of the square. It consists of 8 elements: the 
identity {/}; reflections through the channel axes {a x ,a y }; reflections through the diagonal 
axes {o"i, CT2}; rotations by angle ir/2, -R-7172}; and, rotation by angle ir, {R n }- These 

five sets of elements comprise the five conjugacy classes. It follows that there are five irreps, 
four are one dimensional and one is two dimensional. The character table is given in Table 
1. 

We set out to calculate the five heat kernels corresponding to the five classes. The 
integral corresponding to the identity is 
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Z(I; (3) = — ^ J dxdyd Px dp y e- pn , (2.1) 

where H = (p 2 + p 2 )/2 + x 2 y 2 is the classical Hamiltonian (and the Wigner transform of the 
quantum Hamiltonian). This integral is extensively discussed in Ref. where it is shown 
that it has a logarithmic divergence. We return to this point below. According to Eqs. ( |1.6|) 
and (|1.7| ), the integral corresponding to a y is given by 

1 f 

Z(a y ;(3) = 77^p J dxdydp x dp y e~ ,3H TTh5(y)5(p y ) 

' dx. (2.2) 



V 87rf3h 2 

Since the x integral runs from — oo to oo, this integral diverges even more violently than 
( p.l| ). We will also return to consider this more carefully below. The remaining three class 
heat kernels are well behaved. For the reflection through the diagonal axis ai, we change 
variables to £ = (x + y)/ \[2 and 77 = [x — y)j a/2 so that 

Z(ax;f3) = J d£d'qdp ( .dp ri e- l3H 'nh5('q)5(p ri ) 

- T M^_ (23 ) 

The kernels corresponding to rotations by tt/2 and tt are trivial since all integrals are done 
by delta functions leaving 

Z(R v/2 ;/3) = ± Z(R K ;P) = ± (2.4) 

We now go back and analyse in greater detail the first two integrals. The first was studied by 
Tomsovic [|TU[ but for completeness and consistency of notation, we review the calculation. 

Deep in one of the channels, x ^> 1 for example, the approximation using the Wigner 
transforms breaks down. This can be understood as follows. This approximation assumes 
that for short times one is free to ignore the dynamics so that the calculation involves only 
the local value of the Hamiltonian. Usually this is not problematic, however one finds here 
that the channel effects violate this assumption. This is because in one of the channels, 
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x ^> 1 for example, we can treat the problem adiabatically so that in the y direction there 
is harmonic motion with a frequency uj x = y/2x. This frequency becomes arbitrarily large 
in the channel and there is no time scale over which the dynamics can be ignored. We 
overcome this problem by using an alternate representation of the heat kernels based on the 
approximate separation of the problem as introduced in Ref. JIIJ and which we discuss below. 
This is a complementary representation which is valid deep in the channels but fails near the 
origin. To proceed, we assume that there is a domain of x where both representations are 
valid. Let Q be a value of x in this domain. The condition for the adiabatic representation 
to be valid is that we be deep in one of the channels, in terms of dimensionless quantities this 
is (3 l / A Q ^> 1. The condition for the Wigner function representation to be valid is (3fiQ <C 1. 
These conditions are compatible if P 3 ^ 4 h < 1. (If we determine all quantities in units of 
energy [E], then Q has units of [£" 1//4 ], P has units of [E^ 1 ] and h has units of [i? 3 / 4 ] so that 
all the conditions mentioned above are in terms of dimensionless combinations.) We will use 
the Wigner representation in the square \x\ < Q, \y\ < Q and the adiabatic representation 
elsewhere. We then replace ( |2.1|) by 



(We have introduced the subscript to denote that this is the contribution from the central 
region around the origin.) This integral can be done by the change of variables u = \ffixy 
and v = x, so that the integrand is proportional to exp(— u 2 )/v. Doing the v integral first 
and using (3 l / 4 Q 3> 1, one finds 



where 7 = 0.5772... is Euler's constant. Similarly, for the reflection operator a y , integration 



To do the integrals in the channels we assume a local separation of the Hamiltonian 




(2.5) 




(2.6) 



of ( |2.2| ) between the limits — Q and Q leads to 




(2.7) 
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into a free particle in the x direction and a harmonic oscillator in the y direction, with a 
frequency which depends parametrically on x 

K = \v\ + ^V 2 . (2.8) 

Henceforth we will use small letters to denote objects related to the local Hamiltonian h x . 
It has eigenenergies e n = (n + l/2)u x h and eigenstates \<p n ) which depend parametrically on 
x. All the symmetry information to do with the channel calculation is encoded in these local 
eigenenergies and eigenstates. In particular, we are interested in the subgroup of C± v which 
leaves x invariant and so maps the local eigenstates onto one another. This subgroup is just 
the parity group with group elements {I,a y }. This group has a trivial character table but 
we include it for completeness as Table. 2. For fixed x, we proceed in analogy to (|1.9|) by 
defining heat kernels based on the local eigenvalues and corresponding to these two group 
operations, 

z x (g,(3) = tr(u\g)e^ h ) 

= J2Vn(g)e-^ (2.9) 

n 

where g is either the identity or the reflection element. The trace operator "tr" denotes the 
local integral over the y degree of freedom and can be found by summing over the index 
n. It is clear that the operator W(g) is unity when g = I and changes the sign of the odd 
states when g = <r y , so that r] n (I) = 1 and r] n (a y ) = (—1)". 

To evaluate the full trace, we note that the integrals in p y and y have already been done 
implicitly in ( |2.9| ) so we only need to do the x and p x integrals. Since this is only one 
dimensional, the prefactor of the integral has only one power of 2nh and we conclude 

Z c (g; P) = f 9 ^ 1^ d Px e~^ 2 J q dxz x (g, ft 

where we have defined the factor £ = exp(— phQ/^/2). (We include a subscript c to denote 
that this is the channel contribution.) We have also introduced a factor f g which represents 
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the number of channels which map to themselves under the action of the group element g. 
When working with the identity element, all the channels map onto themselves and // = 4; 
when working with the reflection operator the two such channels along the x axis map onto 
themselves and f ay = 2. We now make use of the series identities 

Y(-l) n - = arctanf (2.11) 

V 2n + l 5 1 ; 

and the fact that (3%Q <C 1 to conclude 



We add the results from inside the square ( |2.6| ) and ( |2.7| ) to the channel results (|2.12| ) to get 



z ^={jk- (2 ' 13) 

Note that the Q dependence has cancelled from both results leaving a finite answer. (This 
prescription actually overcounts some regions of phase space but the errors so introduced 
are exponentially small in /5 1 / 4 Q.) 

We have now calculated the five class heat kernels which we need. All that remains is to 
compute their inverse Laplace transforms. In fact, we will not be interested in the densities 
p(g; E) themselves but rather in their integrals N(g; E) which are given by 

N(g;E) = C-i(^±). (2.14) 



ft 



The inverse Laplace transforms are 

2 



N(I;E) = -V (4 logy + 4 7 + 14 log 2 -8) 
N(a y ,E) = 2 -y 2 
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N(R n/2 ;E) = X - 
N(R W ;E) = ~. (2.15) 

We have defined the dimensionless scaled energy y = E 3 ^ 4 /h, which is a semiclassically 
large quantity. If we explicitly include the mass m in the kinetic energy of the Hamiltonian 
and a parameter a in front of the potential energy then Eq. (|2.15 ) still applies but with 
y = (m 1 / 2 E 3 / 4 ) / (a l l A h) . We further note that the inverse Laplace transforms imply that 
all the functions are zero for negative energies. The first of these relations is the average 
integrated density of states summed over all irreps and was already found by Tomsovic |TU 



To construct the integrated densities of states for each of the five irreps, we use Eq. ( |1.10| ) 
with the symbols Z replaced by the symbols N. It should also be mentioned that these are 
just the leading order results in an asymptotic semiclassical expansion. The terms in this 
series will eventually diverge in a manner controlled by the shortest periodic orbit ||18|| . 

For typical two dimensional potentials with finite phase space volumes, the term N(I; E) 
scales as 1/Ti 2 . The prefactor of that term in ( [2.1 5|) has this scaling but there is a further 
logarithmic dependence on h which causes it to grow somewhat faster. This logarithmic 



factor arises from the fact that the integral in ( |2.5|) diverges logarithmically with Q. One 
must be careful in discussing "orders" when expressions involve logarithms of large quantities 
and for practical purposes, the non-logarithmic term 47 + 14 log 2 — 8 represents an essential 
correction, as discussed in Ref. [Kj. Based on Eq. ( |1.7| ), we expect terms involving reflection 
operators to be weaker by a relative power of h and therefore to scale as 1/h. This is not 
true for N(a y ; E) which is amplified by a factor of 1/h so that it is of the same order as 
the non-logarithmic term in N(I;E). The fact that it has been amplified by a full power 
of 1/h can be traced to the fact that the integral (f2.2|) diverges linearly with Q. Therefore, 
rather than being a relatively weak correction, this reflection operator is almost leading 
order in its effect. In particular, the approximate relation that the fraction of states in 
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irrep R is approximately d 2 R /\G\ fails in general, since it comes from considering just the 
identity operator. (However it is valid for the E irrep which is independent of that reflection 
class.) A similar behaviour is also apparent in the related problem of the hyperbola billiard 
fT^Pfl- The other reflection class function N{a\\ E) does scale as 1/h as we expect for normal 
reflection operations. The two rotation classes also behave normally |12|| , being constants 
independent of Ti. 

III. NUMERICAL COMPARISON OF TWO DIMENSIONAL RESULTS 

We have numerically diagonalised the quantum Hamiltonian and found the first few 
hundred eigenvalues of the problem. We used appropriately symmetrised bases involving 
harmonic oscillator wave functions in the x and y directions to separately find the eigenvalues 
belonging to each irrep. All results are for bases of 200 oscillators in each direction. To make 
the comparison more explicit, we convolved the numerically obtained density of states by a 
Gaussian of width w, 

PR(E) = ^=^£e X p (- iE :f )2 \ (3-1) 

The integrated density of states is then obtained by replacing the sharp steps at the quantum 
eigenvalues by the corresponding error functions. For large w, this convolution washes out 
all oscillations leaving just the average behaviour. 

In Fig. Q we show the results for all five irreps with a smoothing width w = 3. The solid 
curves are the numerics and the dashed curves are the analytical forms. The first thing which 
is apparent is that there is a great distinction between the A\ and B\ states compared to the 
A2 and B2 states, resulting from the large contribution of N(a y ; E). Between each of these 
pairs there is a much smaller splitting due to N(ai]E). The deviations between the solid 
and dashed curves are completely numerical in origin and arise from the finite basis used in 
determining the quantum eigenvalues. Due to the channels, the eigenvalues converges very 
slowly with increased basis size. It is interesting to note that the irreps which are odd with 
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respect to reflections through the channels are better converged. Being odd, they are less 
sensitive to the effects of the channels and are therefore less error prone. Nevertheless, their 
error is still dominated by channel effects as we will demonstrate. The other three irreps are 
not odd with respect to both channels (A± and B\ are even with respect to both channels 
and the E states can be chosen as even with respect to one and odd with respect to the 
other.) All three of them fail at approximately the same energy of E ~ 18. The number 
of accurate eigenvalues is roughly 35 for A\ and B\ and 45 for E (recall that E is doubly 
degenerate so the number of independent eigenvalues obtained is half the number of states 
plotted.) This is rather dismal considering the 40,000 oscillator states used. The irreps A 2 
and B 2 are accurate up to energies near E « 60 representing roughly 115 states each. 

It is also interesting to numerically isolate the contributions from the various classes and 
compare them to ( |2.15|) directly as done in Ref. jl!|. This is a simple exercise since the 



entries in the character table are components of a unitary matrix which is readily inverted. 
The result is 



' N(I;E) N 


(l 


1 


1 


1 


1^ 


' N A1 (E)^ 


N(a y , E) 


1 


-1 


1 


-1 





N A2 (E) 


Niat-E) 


= 1 


1 


-1 


-1 





N m (E) 


N{Rk/ 2 ]E) 


1 


-1 


-1 


1 





N B 2(E) 


v N(R n ;E) / 




1 


1 


1 




v Ne(E) j 



(3.2) 



This can be written compactly as 



N{g-E)=Y j m{g)N R {E) 

R 



(3.3) 

where the factors i]n(g) are defined in ( |3.2j ) and can be thought of as the inverse of the 
group characters. In Fig. we plot N(I; E), N(<j y ; E) and N(R V ; E) from the theory and 
with the numerical eigenvalues combined according to ( |3.2| ). As mentioned, the first is just 
the total number of states. The third is shown in its own panel since its value is of a very 
different scale than the other two. They all fail around E ~ 18 which is consistent with the 
previous figure. N(R n ; E) depends on very fine cancellations and is more sensitive to small 
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errors so it is consistent that it produces noticeable deviations at a slightly smaller energy 
than the other two. Eq. ( |2.15| ) predicts a flat line for N(R n ;E), the structure at smaller E 
comes from the convolution ( |3.1| ) which is applied to the analytical forms as well as to the 
numerical data. 

The other two conjugacy classes behave very differently. We plot these results in Fig. |4j. 
The upper panel shows N(ai;E) and the lower panel shows N(R n / 2 ; E). For the lower 
panel, we choose two different smoothing widths, the relevance of which we discuss below. 
For now, consider the comparison between the smooth solid curve and the dashed curve in 
each case. The results are now accurate up to energies of E « 800 or more than forty times 
the range observed in the previous figure. This indicates that the numerics are, in some 
sense, better than a quick study of Fig. || would indicate. Although the various irreps are 
individually error prone even at relatively modest energies, these errors are very correlated 
so that appropriate combinations cause them to cancel. In fact, this is apparent in Fig. |2| 
since the pairs A\ and B\ and also A 2 and B 2 deviate from their expected behaviour in 
very correlated manners. From Q3.2| ) we see that both N(ai,E) and N{R v / 2l E) involve 
the differences Nai(E) — N B i(E) and Na2(E) — N B2 (E) and the systematic effects cancel 
for these two classes. Since these functions agree with the numerics up to E ~ 800, it is 
reasonable to associate all the problems in the numerics with N(I;E) and N(a y ;E), i.e. 
with the channels. This is obviously true for the irreps Ai, B x and E, however it is also true 
for the odd irres A 2 and B 2 . Their staircase functions fail at E m 60 which is better than 
the other irreps but still very much smaller than the classes N(ax \ E) and N(R n / 2 ; E). 

We now briefly discuss the oscillatory structure visible in the bottom panel of Fig. 01. 



This type of structure was also observed in Ref. |12j where it was explained in terms of 
fractions of periodic orbits ||20|| . In this example, the structure arises from the square-like 
periodic orbit shown in Fig. |T[ After completing, one quarter of a cycle, the trajectory is 
related to its initial point by a rotation of angle n/2. This quarter-orbit then contributes 
an oscillatory contribution to the function N(R n / 2 ; E). This is a scaling system whose 
classical mechanics is independent of energy, after appropriate scalings. In particular, the 
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period of an orbit scales as T a E~ x l^ which explains the growing wavelength with energy. 
Additionally, the smoothing suppresses the oscillatory contribution by a factor proportional 
to exp (—w 2 T 2 /2) which explains why the amplitude of oscillation increases with energy. 
At the highest end of the energy range, one sees the contributions of higher repetitions 
- for example three quarters of the square orbit will also contribute to N{R 7T /2' ) E). The 
function N(R ai ;E) receives contributions from fractional orbits which map to themselves 
under reflection through the diagonal. Examples of this include the diagonal orbit after a 
half period and after a full period. Such structure is visible at the upper end of the energy 
range but is less apparent than in the bottom panel because of the different vertical scale. 

Similar structure exists for the other classes as well but is not visible due to the short 
energy range available. N(R n ; E) receives contributions from one half the diagonal orbit and 
one half the square orbit. N(a y ; E) receives a strong contribution from the almost periodic 
family of orbits corresponding to the adiabatic oscillation deep in the channels (actually, 
from the fractional periodic family which has one half the period.) This is a non-standard 
contribution due to the intermittency, such effects are discussed in Refs. [|]-[7j . The function 
N(I; E) receives contributions from all the complete orbits but not from any fraction of 
them. The periodic orbit theory of this system has been discussed in detail in Ref. |5j and 
the references therein, so we forego a more detailed discussion. 

IV. THE THREE DIMENSIONAL GENERALISATION 

In this section we discuss the three dimensional potential V = x 2 y 2 + y 2 z 2 + z 2 x 2 . This is 
the potential which actually appears in the zero dimensional limit of the SU (2) Yang-Mills 
equations. The symmetry group is that of the octahedral group in which we allow spatial 
inversions — the extended octahedral group. In Fig. || we show a three dimensional constant 
energy contour of the potential and also an octahedron whose vertices are aligned along the 
channel directions. In total there are 48 group elements organised into 10 conjugacy classes. 
This group is the direct product of the inversionless octahedral group and the inversion 
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parity group. The first of these is composed of 24 group elements organised into 5 classes 
PH and we start by enumerating these. First, there is the identity /, which is in a class 
by itself. There is a class of six elements involving rotations by ±7r/2 about any of the 
three axes, such as R x ,n/2- Similarly, there is a class of three elements involving rotations 
by 7r about these axes, for example R Xj7T . There is a class of 8 elements involving rotation 
by ±27r/3 about any of the face-face axes, such as R a ,2-K/z- Finally, there is a class of six 
elements involving rotations by n about any of the the six edge-edge axes, such as Ri^. We 
refer to these classes as G\ to C5 respectively. This group has 5 irreps and the character 
table is the top left quarter of Table 3. 

To construct the full group, we multiply representative members of each class by the the 
inversion operation S = o x a y a z . The effect of this is to map the identity to the inversion 
element E and to map each rotation into either a single reflection or into a rotation times 
a reflection. This induces five additional classes. The element £ is in a class by itself. 
Composition of the second class with £ gives a class of six elements which are rotations 
by ±7r/2 through an axis times reflection through that axis, such as R x ^/2* J x- Composition 
of the third class with £ gives the reflection elements about the three planes, such as a x . 
The fourth class becomes a product of a rotation about a face-face axis times a reflection 
through the perpendicular plane, such as Ra^/z^a- Finally, the fifth class becomes reflec- 
tions through planes defined by the edges and vertices. An example is the plane defined by 
the point 1 together with the vertices at positive and minus z. We call reflections through 
this plane o\. We denote these five additional classes C[ to C 5 respectively. The addition of 
these classes doubles the number of irreps and the full character table is shown in Table 3. 

We proceed by analogy with the two dimensional problem. There we found that to 
analyse the contribution of a single channel, it was necessary to consider the subgroup 
which mapped that channel onto itself — in that case it was the parity group. We do the 
same here. The eight group elements which map the channel x ^> 1 (for example) onto 
itself are /, a VtZ , 02,3, -R x ,±7r/2 and R x ^ and these belong to classes C\, C 3 , C 5 , C2 and C3 
respectively. (02,3 are defined in analogy to a%\ they are reflections through the two planes 
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defined by the vertices at plus and minus x and the midpoints of the two edges connecting 
the z vertex to the positive and negative y vertices.) We can expect the integrals associated 
with these elements to be problematic and to possibly require the adiabatic matching used in 
section III. Together these eight elements comprise the subgroup C± v which is, of course, the 
group we studied in the two dimensional problem. This will prove useful in the subsequent 
analysis. 

We start by studying the five classes which do not require an adiabatic analysis. The 
class C[ involves three orthogonal reflections while the classes C 2 and C4 involve rotations 
and perpendicular reflections. Their Wigner transforms are given by ( |1.8|) and are trivial to 
integrate since they involve delta functions of all the quantities. Their constributions are 
1/8, 1/4 and 1/6 respectively. The class C4 involves rotations through the face axes. For 
rotation by 2ir/3 through the point a, we define a change of variables 

Z=-^(x-y + z) v = ^=(2x + y-z) ( = _( y + z ), (4.1) 

so that the potential along the £ axis is V = £ 4 /3. We then use the third equation of ( |1.7|) 
with this choice of variables to find 

(4 - 2) 

For rotation by tt through the point 1, we define a change of variables 

i = -^{x + y) v = ^(x-y) C = z (4.3) 

so that the potential along the £ axis is V = £ 4 /4. We then find 

We now consider the more interesting classes which map at least one channel onto itself. 
We earlier suggested that the integrals corresponding to them might be problematic. In 
fact, this is true for all of them except the identity whose integral converges without such 
an analysis. Therefore, we do it first, 
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Z(I\ (3) = , 2 J^, 3 J dxdydzdp x dp y dp z e m 



1' 



(4.5) 



(The p integrals are done trivially and the spatial integrals can be done by a change to 
cylindrical coordinates.) The convergence of this integral is due to the fact that deep in one 
of the channels, the energetically accessible area pinches off as 1/x 2 , which is integrable. 
The analogous integral in two dimensions pinches off as 1/x and is not integrable. The 
remaining four classes follow from using (|1.6|), (|1.7|) and ( |1.9|) inside a cube \x\ < Q , \y\ < Q 
and \z\ < Q. 

For reflection in z, which is a member of the C ; 3 class, we use Eq. ( |1.7| ) and so arrive at 
the following integral, 

I 1 poo rQ 

M^P) = 7:77^2 / dp x dp y e-^ + ^ 2 / dxdye-^ 2 . (4.6) 



2 (27T^) 2 J-00 " ry J-Q 

Other than the factor of one half, this is the same integral we evaluated to get the total 
density of states in the two dimensional problem. The result is given by fl2.6| ) and so we 
conclude 



Z °(^ P) = ^ ^4 (^g Q + log/? 1/4 + J + \ log 2) . (4.7) 

Reflection in 03, which is a member of the C' 5 class, requires a more complicated calculation. 
We define a change of coordinates so that r\ = (z + y)/ \pl and ( = (z — y) / \/2 and then use 
Eq. ( |1.7|) with the delta functions acting on ( and p^ so that the integral to be evaluated is 

Z (a 3] (3) = — -y / dx ^e-^V + W4). (4g) 
rcph Jo Jo 

We have done the trivial momentum integrals and have noted that by its definition, r] has 
a different integration range than x. This integral can be done in a manner analogous to 
Q2.5| ), we define integration variables u = xrj and v = rj. Doing the v integration first and 
using /3 1//4 Q ^> 1 one arrives at 



Zo{a3] ® = V A^Pf? ( l0gQ + l0g/?1/4 + \ + I bg 2 } • (4 ' 9) 



Rotation by 7r/2 about the x axis is a member of the C 2 class and implies delta functions 
in the other two variables so that the integral to be done is 



Z (R 7r/2 ;P) = — dp x e-^' 2 dx 
Ann J-00 J-o 



~Q. (4.10) 



V 2nph 

Rotation by it about the x axis, which is a member of the C3 class, involves an integral 
which is identical except for a factor of two from the sin 2 (0/2) factor in (|1.7|) . Therefore 



z ° ( ^> = (411) 



V. CHANNEL CALCULATIONS IN THREE DIMENSIONS 

In this section we evaluate the contribution of the channels in three dimensions. As 
discussed before, this is is only necessary for some of the group elements. In analogy with 
Q2.8p we define a local two-dimensional Hamiltonian as 

K = l(pl+ P 2 z ) + ^(y 2 + z*) + y^ } (5.1) 

where again to x = \/2x and x is assumed large. Deep in the channel, the final term can 
be thought of as a small perturbation which has virtually no effect on the eigenenergies. If 
that term were completely absent, the local Hamiltonian would have an SU(2) symmetry 
corresponding to a two-dimensional harmonic oscillator. The eigenvalues of the Hamiltonian 
would then be e n = (n + l)hu x , each with a degeneracy of (n + 1). The degenerate states 
can be labelled by the rotational quantum number m which runs from — n to n in even 
increments. The perturbation y 2 z 2 will not affect the energies in a significant manner but 
will act to break up the degenerate collections of states into specific irreps of C± v as follows. 
All states with odd m correspond to the E irrep. The m = states are all A\. For m 
non-zero and divisible by 4, the states are either A\ or B 2 (corresponding to cos(m#) and 
sin(m#) respectively). Otherwise, if m is even but not divisible by 4, the states are either 
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A2 or B\ (corresponding to sm{md) and cos(m8) respectively.) We then define local heat 
kernels corresponding to the five irreps by adding the contributions of all values of n with 
the appropriate degeneracy factor for each irrep so that 

~n + 4" 



zbM 
z Bl ((3) = z A2 ((3) 



E 

n=even 

E 

n=even 

E 

n=cvcn 



-j3hu) x (n+l) 



4 



-f3hu) x (n+l) 



n=odd 



(5.2) 



where [x] is the largest integer less than or equal to x. We will refer to these relations 
collectively as 



-/37kj x (n+l) 



z r(P) = E c R( n ) e 

n=0 

where cr(h) are the degeneracy factors defined in 

To evaluate the traces, we integrate over the remaining x dependence 

Z R ((3) = — ^c R (n) / dp x e~^ 2 / dxe~^ n+1 > 
LtxTi „ J-00 Jo 



(5.3) 



n 



n+l 



(5.4) 



where we have defined ^ = exp (-V2/3HQ) w 1 - V2phQ. (Note that this is different by 
a factor of two from the analogous variable in two dimensions.) All of this discussion is 
in terms of the local irreps; what we really want, however, are the local class heat kernels. 
These we can get by appropriate combinations of the irreps as in ( |3.3| ) to arrive at the class 
sums 



%;/3) = 



71 



n + l 



(5.5) 



For the moment we omit the prefactor of ( |5.4| ), this will be reintroduced later. The de- 
generacy factor c(g, n) corresponding to a group element g is found by adding together the 
degeneracy factors cr{ti) with the appropriate weightings as given by 



i.e. 
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c(g,n) = J2r] R (g)c R (n). (5.6) 

R 

We start with the identity element. Earlier it was argued that we do not need a channel 
calculation since the central integration converges. However, it is of interest to see how this 
is also apparent in the channel calculation. Comparing ( |3.2|) and ( |5.2|) it is apparent that 
c(I, n) = n + 1 so that 

OO £ 

S(i; (3) = £ T +1 = t^t 

n=0 1 S 

~ y/2P%Q' ^ 



We now need to reinsert the prefactor of (|5.4j ) and must also include an integral factor 



representing the number of channels left invariant by the corresponding element f g , as in 
two dimensions. We trivially have fj = 6 so that the channel result for the identity element 
is 

z ^-1T,WWq- (5 ' 8) 

We now compare this result to ( |4.5| ), the present contribution is very much smaller if /3 1 / 4 Q S> 
1 which is precisely the limit we are considering. Therefore, we again observe that no channel 
calculation is necessary for the identity element. 

We next consider the reflection element o z . It is in the same class as a y so comparing 
( P-2| ) and (|5.2p we conclude c(a z ,n) = 1 when n is even and when n is odd. The sum 
which must be done is 

tn+l oo c2m+l 

S(a z ;(3)= Y i 

n=even n < 1 m=o Z '" L ' 1 



where we have used ( 2.11| ) and the approximation immediately below (|5.4j ). Note that 



= 4 since a z leaves four channels invariant, so that 



21 



Recalling now the corresponding result for the central region (|4.7| ), we conclude that for the 
class 



z(c ^Wi^H^) +7+4ios2 )- (5 ' n) 

This is independent of Q as we expect. 

The equality of zb 1 and za 2 in (|4.7|) implies that S{a^[3) = S{cr z ;(3) (since they both 



equal za 1 — zb 2 from ( |3.2| ).) The only difference is the subsequent calculation is that f as = 2 
so that the channel calculation is one half of that for a z ( |5.10| ). We combine this result with 
the result from the central region ( f4.9|) to determine 



z(C ^W«H^) +7+81 °4 (5 - i2) 

For rotations by 7r/2 about the x axis we note that c(i? vr / 2 , n) = (— l) n / 2 for n even and 
is for n odd so that 

4 



n+l 00 £2m+l 

m 



cn+l 00 

S(R n/2 ;P)= £ ("ir /2 f— =E("1 



ir. OVOU " ' n>- H ^ m + 



arctan£, (5.13) 



where we have again used ( j2.11| ). We now note that arctan£ « 7r/4 — (3hQ/\/2 and also 



that only two channels are left invariant implying /r. 2 = 2 so that 



-./i^-V^e- (5 - 14) 



We now combine this with the calculation from the central region ( 4.10|) to arrive at 



z ^ )= iwff- (5 - 15) 

The final class to be analysed is C3 of which rotations about the x axis by 7r is a 
representative member. We now have c(R n ,n) = (—1)™ so that the relevant sum is 

OO £ 

S(R 7r ;(3) = J2(-m n+1 = T^ 7 

-sHM- (516) 
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As for the previous rotation class, we have fj^ = 2 so that the result of the channel calcu- 
lation is 



Combining this with the calculation from the central region ( |4.11| ) we conclude 



z(c ^ )= i^w (5 - 18) 

The final analysis we will do is to find the inverse Laplace transform of the various 
relations and thereby express them in the energy domain. The ten results as a function of 
(3 are scattered over the previous two sections. As in the two dimensions, we go directly to 
the integrated densities of states by use of ( |2.14j ). The result is 



1&T 2 (±) 
N(C i; E) = ^My 3 



N(C 2 ;E) = \y 2 



N(C 3 ;E) = ^-y 2 



45%/^ 
1 

3' 
2 

3V 



N{C 4 ;E) 



r 2 (|; 



27^3tt 3 



1 

8 
1 
4 

N(C' 3 ;E) = -^/ 2 (41og2/ + 4 7 + 101og2-8) 

D7T 



N(C[;E) 
N(C' 2 ;E) 



N{C' A -E) 



1 
6 
1 



N(C' 5 ;E) = -V(41og<, + 4 7 + 141og2-8), (5.19) 

where again we use the semiclassically large quantity y = E 3 '^ /%. For comparison, we 
remark that for generic potentials, use of ( |1.7j ), would imply that the first term scales as y 3 , 
the following four as y, the set {C[, C' 2 , C'^} as y°, and the set {C' 3 ,C 5 } as y 2 . 

23 



The leading order behaviour, as given by the first expression, scales generically with %. 
There are no other terms which are competitive with it so the relation that the fraction of 
states in irrep R is approximately d 2 R /\G\ is valid. As discussed in the text the reflection 
classes C' 3 and C' 5 are amplified somewhat, having an additional logarithmic dependence on 
fi in addition to the 1/h 2 prefactor. This is in analogy to the total density of states of the 
two dimensional problem. In fact, the class C' b is, within a factor of four, the same as the 
total density of states in two dimensions. Two of the rotation classes are amplified by 1/h so 
that they scale as 1/h 2 . This makes them competitive with the reflection classes (since, as 
argued in the two dimensional problem, the logarithmic term is a rather weak amplification). 
This is analogous to the behaviour of one of the reflection operators in the two dimensional 
case. 

In Fig. |6| we show the integrated densities of states found from using the results of (|5.19|) 
combined according to the characters of Table III. (It should be remarked that this is not 
entirely consistent since the leading order terms have semiclassical corrections which are 
almost certainly of the same order or larger than the smallest terms we are considering. 
However, the point of this paper is not a systematic semiclassical expansion but rather a 
study of the symmetry effects.) The structure now looks more typical; irreps of the same 
dimensionality have roughly similar numbers of states with slight differences arising from 
the contributions of the other group elements. In particular, the largest four curves are the 
four three dimensional irreps and the differences among them arise from the terms of order 
y 2 log y and y 2 ; the largest of these curves belongs to r' 5 . The middle two curves belong 
to the two dimensional irreps and the smallest four curves belong to the one dimensional 
irreps. The largest of these is the trivial irrep Fx; this is reasonable since it receives positive 
contributions from all the classes. 

In Fig. ^ we show the same data but on a smaller energy scale. At the right edge 
of the figure (E = 35), the curves are ordered the same as in Fig. || (i.e. their asymptotic 
ordering). However, it is clear that there is a lot of crossing of these curves at lower energies. 
This is because for moderate energies the contribution corresponding to identity in (|5.19|) 
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does not dominate the others. Additionally, in calculating the functions for each irrep via 
( P-.10| ) (with the symbol Z replaced by N), we must sum over all the group elements and 
so the contribution of any given class is amplified by the number of elements in that class. 
The identity class only has one element but the classes which contribute to next order, 
{C2,C 3 ,C' 3 ,C' 5 }, Yi&ve six, three, three and six elements respectively. As mentioned, it is 
difficult to calculate many accurate eigenvalues when a potential has channels and this is 
especially true in three dimensions. Therefore, the non-asymptotic behaviour in Fig. [7| is 
relevant to any numerical study since the results will probably all be in that energy domain. 



CONCLUSION 



We have shown that the symmetry reduction of the Thomas-Fermi density of states 



discussed in Ref. is easily generalised to more perverse systems where the Wigner rep- 
resentation fails. In two dimensions, the symmetry decomposition introduces essentially 
leading order contributions to the densities of states of the one dimensional irreps. The 
results were verified numerically and seen to work well. However, the problem studied is 
numerically very difficult and only a handful of states of each irrep are reliably calculated. 
Nevertheless, certain combinations of the densities of states are found to be accurate to very 
high energies even though the density of states of each individual irrep is not. This effect is 
noticeable only by studying the class functions derived here and would not otherwise have 
been apparent, thus underlining the importance of symmetry decompositions. 

In three dimensions, we find that the symmetry decomposition does not introduce terms 
which are essentially leading order. However, there are still interesting effects; two of the 
reflection classes have a logarithmic dependence on fi beyond what one might have expected 
and two of the rotation classes have an additional power of 1/h thus making them of es- 
sentially the same order as the reflection elements. Furthermore, we observed that even in 
this case one must consider rather high energies before the ordering of the functions Nr(E) 
achieves its final form. This is in spite of the fact that the leading behaviour is not affected 
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by the decomposition. Rather it arises from the fact that the classes which contribute at 
next to leading order have several group elements and their contributions are correspond- 
ingly amplified. This is an effect which we can expect to become even more important in 
higher dimensions if we consider potentials of the form V({xi}) = J2iJ2j>i^x^. In higher 
dimensions, more and more of the terms will behave with the normal ft dependence. The 
only terms with anamalous dependences are those for which one would initially expect a 
dependence of 1/ft 2 or 1/ft. If the corresponding group element leaves at least one channel 
invariant, they will be amplified by factors of log (1 /ft) and 1/ft respectively. 
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FIGURES 

FIG. 1. The configuration space of the x 2 y 2 potential. The light solid lines are constant energy 
contours at E = 0.1,1,3,5,--- and the four dashed lines indicate the axes through which the 
system has a reflection symmetry. The two heavy curves show the two shortest periodic orbits in 
the system calculated at E = 1. 

FIG. 2. The solid curves indicate the smoothed density of states for each of the five irreps as 
found numerically. The dashed curves are the corresponding analytical forms derived in this paper. 

FIG. 3. Top: The upper pair of curves indicates the function N(I; E) which is the theory for 
the total density of states. The lower pair shows N(a y ; E). In each case the solid curve comes from 
the numerics and the dashed curve is the analytical form. Bottom: the same for N(R W ; E). 

FIG. 4. The same as the previous figure except that the upper panel indicates N(a\;E) and 
the lower panel is N(R 7r / 2 ; E). The solid oscillating curve has with a smoothing width of 3. The 
other solid curve and analytical dashed curve have smoothing widths of 30. 

FIG. 5. Left: An equal energy contour of the three dimensional potential 

V = x 2 y 2 + y 2 z 2 + z 2 x 2 showing the six channels along the three axes. Right: An octahedron 
with the relevant points labelled for the description of the group elements. 

FIG. 6. The average densities of states for each of the ten irreps of the V = x 2 y 2 + y 2 z 2 + z 2 x 2 
potential. From greatest to smallest the curves describe the irreps r' 5 , T4, T^, F$, T3, T' 3 , T±, T' 2 , 

r 2 and r;. 

FIG. 7. The same as Fig. 6 but on a smaller energy scale to show the curves crossing at small 
energies. At E = 35, the order of the curves is the same as that described in the previous figure 
caption. 
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TABLE II. Character table of the parity group. 
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TABLE III. Character table of the extended octahedral group. The number in brackets at the 
top of each column indicates the number of group elements which belong to that class. Represen- 
tative members of the various classes are described in the text. 
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